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SUMMARY 


This  volume  presents  numerical  calculations  of  ionospheric 
electron  density  perturbations  and  ground-level  signal  changes  for 
several  electric  field  distributions  produced  by  high-frequency  (HF) 
transmitters.  Plasma  density  perturbations  are  likely  to  occur  in 
regions  of  intense  electric  field  focusing.  These  shifts  in  the 
dcr.aitv  ci.d  di  sti  iijuLion  or  the  plasma  near  the  reflectrix  mav  be 
manifested  in  measurable  ground-level  signal  changes.  We  have  simu¬ 
lated  these  interactions  in  our  calculations. 

Our  analysis  takes  into  account  radio  field  focusing  at  caustics, 
fhe  consequent  j ou le-hea t ing  of  the  surrounding  plasma,  heat  conduc¬ 
tion,  diffusion,  and  recombination  processes--these  being  the  effects 
of  a  powerful  oblique  "modifying"  wave.  We  then  seek  effects  on  a 
secondary  "test"  wave  which  is  propagated  along  the  same  path  as  the 
first.  The  test  wave  could,  in  fact,  be  the  modifying  wave  itself,  in 
v;hich  case  the  calculation  would  be  of  the  "self-effect"  of  a  powerful 
transmitter . 

Our  calculations  indicate  that  ground-level  field- intensity 
changes  of  several  dB  might  be  produced  by  joule-heating  of  the  iono¬ 
sphere  bv  intense  oblique  HF  waves.  Although  small,  these  calculated 
changes  are  similar  to  the  3  dB  change  measured  by  Bochkarev  and 
associates  [1,2.3,‘4],  who  gave  only  sketchy  details  of  their  ex¬ 
perimental  parameters.  Our  results  are  extremely  sensitive  to  the 
model  ionosphere  and  therefore  indicate  that  an  experiment  should 
employ  the  widest  possible  range  of  frequencies  and  propagation 
conditions.  An  effective  power  of  90  dBW  is  far  more  likelv  to 


produce  a  detectable  signal  change  than  85  dBW. 
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PREFACE 


This  volume  applies  theory  developed  in  earlier  work*  [5]  to  the 
problem  of  possible  self-effects  of  powerful  obliquely  incident  radio 
waves.  Using  the  results  from  this  earlier  work,  we  use  field-driven 
changes  in  ionospheric  electron  density  to  calculate  modified  ambient 
refractive  index  profiles  and  the  consequent  effects  on  ground  signal 
intensity  for  sample  transmitters.  This  volume  is  not  self-contained 
and  should  be  used  in  conjunction  with  Vol.  1,  which  presents  the 
underlying  equations  and  defines  the  majority  of  the  symbols. 


’'Field,  E.  C..  R.  M.  Bloom,  and  K.  E.  Heikes.  Ionospheric  Heating 
with  Oblique  Waves.  Vol.  1.  Electron  Density  Perturbations. 
Pacific-Sierra  Research  Corporation.  Report  1864,  September  1988, 
AFGL-TR-88-0336.  ADA216260 
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I.  INTRODUCTION 


Volume  I  of  this  report  1  ’  ]  derived  equations  for  the  electron 
density  perturbations  produced  bv  an  oblique  angle  ionospheric  heater, 
and  presented  some  results  for  idealized  heating  functions  that  al¬ 
lowed  analytic  solutions.  Tliis  volume  presents  numerical  calculations 
of  ionospheric  electron  density  perturbations  and  ground-level  signal 
changes  for  several  spatial  electric  field  distributions  produced  by 
HF  transmitters.  Ve  have  assu.med  antenna  patterns  and  power-gain 
products  based  on  large  shortwave  broadcast  transmitters,  such  as  the 
new  Voice  of  America  facilities  (e.g.,  (6)).  We  employ  simplified 

daytime  and  nighttime  two-laver  ionospheric  profiles. 

It  is  the  nature  of  oblique  propagation  that  electric  fields 
focus  strongly  into  caustics  'which  compensate  for  signal  loss  due  to 
geometric  spreading  (7,81.  The  significant  region  of  intense  focusing 
occurs  along  the  Rav-Bundle.  just  beyond  the  wave-group  reflection 
height  in  the  ionosphere.  l-Jhat  plasma  density  perturbations  may  occur 
as  a  result  of  the  field  are  therefore  likely  to  be  most  pronounced  in 
tlic  cicini'.,,  of  this  It  i.^  loasor.able  to  assume  that  even 

subtle  shifts  in  the  density  and  dimension  of  the  plasma  near  the 
reflectrix  will  be  manifested  in  more  profound  effects  in  ground- 
level  signal  strength--in  particular,  at  the  edge  of  the  geometrical 
shadow.  It  is  this  effect  which  our  calculation';  are  detlr''pd  *"0 
measure . 

Our  analysis  takes  into  account  radio  field  focusing  at  caustics, 
the  consequent  joule-heating  of  the  surrounding  plasma,  heat  conduc¬ 
tion,  diffusion,  and  recombination  processes--these  being  the  effects 
of  a  powerful  oblique  "modifying"  wave.  We  then  seek  effects  on  a 
secondary  "test"  wave,  which  is  propagated  along  the  same  path  as  the 
first.  The  test  wave  could,  in  fact,  be  the  modifying  wave  itself,  in 
which  case  the  calculation  would  be  of  the  "self-effect"  of  a  powerful 
transmitter . 

In  essence,  our  approach  means  recalculating  an  electric  field 
distribution  using  an  index  of  refraction  modified  by  the  intense 
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wave.  Our  calculations  are  grounded  in  what  might  be  called  a  "first- 
order"  approximation  to  nonlinear  wave  propagation.  We  have, 
moreover,  omitted  from  the  analysis  dependence  of  kinetic  and  recom¬ 
bination  coefficients  on  temperature.  The  controlling  hypothesis  is 
that,  should  first-order  calculations  fail  to  produce  notable  heating 
effects,  nonlinear  effects  on  collision  frequency  and  diffusion  would 
probably  not  increase  the  magnitude  of  the  heating  to  a  detectable 
level.  Certain  nonlinear  phenomena,  such  as  plasma  instabilities, 
might  produce  detectable  effects  omitted  from  the  linear  theory. 
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II.  TECHNIQUE  OF  CALCULATION 


Our  cnlrula:;  ions  are  accomplished  through  an  aggregation  of  three 
separate  procedures,  each  producing  a  .set  of  data  utilized  as  input  bv 
the  next  procedure.  The  first,  uhich  we  shall  call  "Caustic"  !9], 
computes  an  electric  field  distribution,  given  an  oblique  transmitter 
and  a  horizontallv  .stiuctified  ionosplieric  profile.  The  CLiustic  cal¬ 
culation  evaluates  the  <‘lectric  field  as  a  superposition  of  quasi¬ 
plane-wave  components,  each  component  being  an  asvmptotic  solution  to 
"he  wave  equation  iti  a  medium  with  specified  '  f  ie  Id- i  ndependent ) 
retractive  index.  The  I'austic  calculation,  unlike  standard  first- 
order  k'KB  asvmptot  ics .  is  uniformlv  valid,  that  is.  both  in  the  region 
of  geometrical  opcic.s  ;ind  across  caustics.  This  method  allows  us  to 
calculate  field  magnit'.ides  throughout  large  regions  of  space  without 
having  to  piece  together  boundary  conditions  and  different  asymptotics 
inside  and  outside  areas  of  strong  focusing.  When  we  calculate  the 
wave  field  resulting  from  refraction  through  a  perturbed,  no  longer 
horizontallv  stratified  ionosphere,  a  less  accurate  method  than  Caus¬ 
tic  must  be  used  to  ostiivaCe  field  intensitv. 

The  second  procedure  is  the  "Diffusion"  calculation,  utilizing 
the  steady-state  solution  Eq .  (19)'  for  electron  density 
perturbation.  The  electric  field  intensity  data  Efx.z)  from  Caustic 
is  converted  to  data  E(x.z’),  where  r.’  designates  a  coordinate  in  the 
direction  of  the  local  magnetic  field;  we  have  been  examining, 
primarily,  east-west  propagation,  and  thus  z'  lies  iii  a  plane  perpen¬ 
dicular  to  the  direction  of  propagation.  It  is  assumed  that  the  only 
significant  plasma  diffusion  will  occur  in  the  z’  direction.  Average 
values  for  the  coefficients  ,  Ly.  etc.  in  Eq .  i  I'M  are  obtained  from 
empirical  models  of  the  height  dependence  of  ambient  electron  recoin- 
l)ination  coefficients  ind  electron  temperature  (IN],  Equation  (I'M  is 
then  evaluated  for  n/Nr).  the  relative  change  in  electrcii  density,  bv 
numerical  integration.  In  an  identical  manner,  the  data  E(x.z')  can 


''All  equation  numbers  lower  than  A6  refer  to  equations  in  Vol.  1. 
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be  used  to  evaluate  the  relative  change  in  electron  temperature  AT/T^ 
by  numerically  integrating  Eq .  (8).  It  is  Eq .  (19)  and  o/Nq,  chough, 
Chat  ultimately  we  are  interested  in  applying. 

The  last  procedure  in  the  calculation  is  the  estimation  of  the 
ground-level  electric  field  intensity  change  (in  the  "test"  wave)  due 
to  ionospheric  perturbation  o/Mq.  As  mentioned  above,  the  "Caustic" 
calculation  cannot  accommodate  ionospheric  profiles  that  have 
gradients  in  other  than  the  vertical  (z)  direction.  This  is  now  the 
case:  we  are  given  to  consider  Che  propagation  problem  through  a 

medium  with  "modified"  complex  polarizability  a  (x.  z)  = 

'7q(z)(1  +  o/Nq) .  It  is  here  that  rav  tracing  and  geometrical  optics 
calculations  serve  best.  A  well-tested  rav  tracing  program  [11]  is 
used,  first  with  refractive  index  N(z),  then  with  N'(x,z).  In  con¬ 
junction  with  the  rav  trace,  we  integrate  certain  quantities  measuring 
Che  rate  of  change  of  ray  coordinate  with  respect  to  infinitesimal 
changes  in  ray  initial  conditions  (i.e.,  launch  and  azimuth).  This 
allows  for  the  continuous  calculation  of  the  relative  convergence  of  a 
small  (infinitesimal)  pencil  of  rays  around  the  main  ray.  This  con¬ 
vergence  (or  divergence)  factor  is  proportional  to  the  electric  field 
intensity  at  a  point.  Conservation  of  energy  along  ray  tubes  jus¬ 
tifies  this  calculation.  The  product  of  the  component  of  the  Povnting 
vector  in  the  direction  of  local  group  velocity  and  an  infinitesimal 
rav  tube  area  is  constant  over  the  rav  path.  As  rav  tube  area  shrinks 
with  focusing,  the  field  intensity  grows  proportionally.  This  cal¬ 
culation  in  fact,  gives  the  slowly  varying  amplitude  factor  of  the 
geometrical  optics  theory.  We  have  followed  the  methods  of  Nickisch 
and  Buckley  [12,13]. 

We  make  this  calculation  on  the  ground  for  both  the  regular  and 
perturbed  ionosphere  and  obtain  the  chanE,e  in  field  intensity. 

Several  points  must  be  .stressed.  First,  ^lie  geometrical  'iptics  field 
intensity  becomes  unbounded  at  the  skip  distance,  both  at  the  interior 
caustic  and  at  the  horizon.  Analysis  shows  that  the  region  of 
validity  is  wherever  the  tract lonal  change  in  adjacent  rav  spacing 
In  one  wavelength  is  much  less  than  one.  The  method  mav  exhibit  a 
shift  in  the  interior  skip  distance,  but  cannot  predict  accurate 
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fields  there,  which  may  change  on  the  order  of  5  to  10  db .  Second,  in 
the  illuminated  region  past  the  interior  skip,  there  are  typically  two 
rays  landing  at  each  point.  This  gives  rise  to  rapidly  varying  inter¬ 
ference  phenomena  in  the  field  amplitude  depending  strictly  on  the 
phase  difference  between  the  high  and  low  rays  at  each  point.  The 
amplitude  of  the  complete  field  lies  in  the  envelope  formed  by  the  sum 
and  difference  of  the  component  ray  amplitudes.  It  may  be  assumed 
that  the  resultant  amplitude  is.  on  the  average,  the  root -mean- square 
of  the  two  components. 

It  should  be  noted  that  a  related,  though  different,  approach  to 
weakly  nonlinear  wave  propagation  can  be  found  in  work  by  Bochkarev 
[ 1 , 2 , 3 ,  <4  ] .  Here,  the  nonlinearity  in  dielectric  constant  is  assumed, 
a  priori,  to  be  proportional  to  the  electric  field  intensity  lEK. 

That  is,  N ' 2  +  a  I E 1 2 ,  where  N  is  the  refractive  index  and  a  is 

the  coefficient  of  nonlinearity.  Except  for  the  region  immediately 
surrounding  the  caustic,  the  field  is  found  bv  classical  geometrical 
optics  with  the  ambient  refractive  index  N.  The  classical  field  is 
used  as  a  boundary  condition  on  a  small  rectangular  domain  around  the 
caustic;  therein,  a  technique  known  as  the  Method  of  Parabolic 
Equations  is  used  to  approximate  the  field  close  to  rays  [10].  It 
would  appear  that  our  method  of  solution  has  the  advantage  of  provid¬ 
ing  insight  into  the  constitution  of  the  nonlinearity  that  Bochkarev 
has  bottled  up  into  the  ad  hoc  constant  a. 

For  a  treatment  analogous  to  ours,  with  a  deeper  analysis  of  HF 
wave  energy  deposition,  as  well  as  estimates  of  nonlinear  absorption 
due  to  plasma  instabilities,  see  Meltz,  Holway,  and  Tomljanovich 
(14).  Their  work  applies  only  to  vertically  incident  waves.  For  an 
adaptation  of  this  class  of  calculations  to  an  automatic,  feedback- 
propelled  computational  scheme,  see  Bernhardt  and  Dincan  [15). 
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III.  EXAMPLES 


We  distinguish  two  calculations;  the  curtain  antenna  case, 
modeled  to  roughly  approximate  large  HE  facilities  such  as  future  VOA 
sites  [6];  and  the  Bochkarev  model,  extracted  from  a  series  of  ar¬ 
ticles  [1,2, 3, 4]  whose  results  we  have  used  (with  some  skepticism)  as 
a  guide.  Both  day  and  nighttime  ionospheres  are  treated  for  the 
curtain  antenna  model;  whereas  for  the  Bochkarev  model  we  use  only  a 
daytime  model  ionosphere--as  close  to  Bochkarev's  as  could  be  ascer¬ 
tained  from  the  published  account.  In  both  cases,  propagation  is 
perpendicular  to  the  magnetic  field  vector;  the  dip  angle  is  60°. 
corresponding  to  the  middie  latitudes. 

The  antenna  patterns,  transmitter  power,  and  ionospheric  profiles 
distinguishing  these  cases  are  shown  in  Figs,  la  and  lb.  Sample  ray 
traces,  one  for  each  profile,  are  shown  in  Figs.  2a  through  2c.  The 
fact  that  the  curtain  antenna  profile  is  denser  at  its  maximum  (F2 
layer  at  300  km)  than  the  Bochkarev  profile  (F2  layer  at  250  km)  means 
that  the  maximum  usable  frequency  (MUF)  for  the  given  launch  angle  for 
the  curtain  antenna  will  be  greater.  Though  the  MUF  for  the  Bochkarev 
model  was  not  given  in  the  published  articles  [3,4].  by  trial  and 
error  we  found  that  15  MHz  was  appropriate  for  the  propagation  path 
that  was  shown  in  the  articles.  We  found  that  any  degree  of  ar¬ 
bitrariness  in  the  choice  of  ionospheric  profile  makes  all  but  the 
grossest  characteristics  of  subsequent  calculations  fluctuate 
substantially.  In  light  of  this  strong  sensitivity  to  model 
parameters,  we  are  justified  in  abstracting  from  our  calculations  onlv 
a  few  general  principles. 

Figures  3a  through  3c  show  the  calculated  electric  field  strength 
contours  for  the  three  cases.  Each  diagram  is  centered  about  the  wave 
reflection  height  and  caustic  region:  this  can  be  seen  by  comparison 
with  the  corresponding  ray  traces.  All  of  the  indicated  cases  produce 
field  strengths  as  large  as  several  tenths  of  a  volt  per  meter.  The 
diffusion  calculation  then  proceeds  with  the  electric  field  data  and 
various  altitude-dependent  parameters  (see  Table  2  in  Vol.  1, 

Sec.  2). 
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a.  Curtain  antenna,  day. 


f  =  17.75  MHz 
21.5  dB  maximum 
1.25  X  10®  watts 
80  dBW  effective  power 


b.  Curtain  antenna,  night. 


f  =  9.75  MHz 
25  dB  maximum 
10®  watts 

85  dBW  effective  power 


c.  Bochkarev  model. 


f  =  15MHz 
25  dB  maximum 
10®  watts 

85  dBW  effective  power 


Figure  la.  Model  antenna  patterns. 
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Figure  lb.  Model  ionospheric  profiles 
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is  of  primary  interest,  as  the  field  E  appears  only  in  the  combination 
E/Ep  (Eq.  [19]).  We  have  used  nominal  values  of  electron  temperature 
Tg  versus  altitude.  To  calculate  6,  the  fractional  loss  of  energy  per 
electron  per  collision,  we  have  used  the  formulations  given  in 
Table  1.  The  term  (E/Ep)^  can  be  regarded  as  proportional  to  the 
energy  in  the  propagating  wave  lost  to  kinetic  heating  of  free 
electrons.  At  night  Ep  is  somewhat  larger  than  in  the  daytime, 
despite  its  proportionality  to  uj;  that  behavior  occurs  because  the 
quantity  6  comes  to  be  dominated  by  collisions  with  atomic  oxvgen  at 
night  (as  opposed  to  elastic  collisions  with  ions  in  the  daytime). 

Table  2  gives  nominal  parameter  values.  In  order  to  perform  the 
calculations  we  interpolated  from  the  values  in  Table  1  to  obtain 
values  near  the  altitude  of  largest  electric  field.  Figures  4a  and  4b 
show  electron  temperature  changes  caused  by  the  electric  fields  in 
Figs.  3a  and  3b.  The  maximum  temperature  changes  are  on  the  order  of 
a  few  percent  for  the  curtain  antenna  daytime  model,  and  a  few  tens  of 
percent  for  the  Bochkarev  model.  The  difference  in  results  fi^r  those 
two  models  indicates  the  extreme  sensitivity  to  ionospheric  parameters 
and  transmitter  characteristics. 

Figures  5a  through  5c  show  the  calculated  relative  electron 
density  perturbations,  n/Ng,  which  are  not  larger  than  a  few  percent, 
and  usually  are  smaller  than  that.  Those  values  can  be  compared  to 
the  estimates  obtained  in  Vol.  1.  Sec.  IV,  Eqs.  (38)  through  (44K  for 
an  idealized  square  heated  region.  Using  Fig.  3c  as  an  example,  the 
most  intensely  heated  spots  extend  over  dimensions  of  several  tens  of 
kilometers,  whereas  (E/Ep)"  is  on  the  order  of  0.30.  From  Eq .  (44)  we 
have  the  estimate: 
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TABLE  1.  Parameters  For  Diffusion  Calculation. 
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TABLE  2.  Ionospheric  Parameters. 
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Fractional  change  in  electron  density,  Bochkarev 
model ,  15  MHz,  day. 
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where  a  is  the  width  of  the  heated  region:  and  Lfj .  L-j- ,  and  R-p  are 
taken  from  Table  2.  Vol.  1.  Thus,  taking  a  =  10  km,  and,  for  ex¬ 
ample,  values  of  Lfj ,  Lp ,  and  Rj  for  250  km  (9.2  km.  25  km,  and  .54), 
we  see  that  this  a  priori  estimate  is  2  percent,  which  is  close  to  our 
numericallv  calculated  range  of  h/Nq  shown  iit  Fig.  5c.  Wliile  the  iso- 
contours  of  field  strength  conform  to  the  downward-curving  pattern  of 
tlie  caustic,  those  of  the  density  perturbation  are  elongated  in 
roughly  a  vertical  direction.  This  is  because  of  our  assumption  that 
diffusion  occurs  onlv  in  the  direction  of  the  magnetic  field  lines 
(see  Fig.  6) . 

Figure  7  shows  the  final  result  of  the  calculation.  This  is  the 
change  in  the  geometrical  optics  field  intensity  at  the  ground  caused 
by  the  heater-induced  ionospheric  perturbation.  Phase  as  a  function 
of  ground  range  is  obtained  by  integrating  the  refractive  index  over 
the  ray  path.  Power  density  as  a  function  of  ground  range  is  calcu¬ 
lated  by  Nickisch's  method  (12).  There  are  generally  two  or  more 
rays,  each  with  its  own  pliase  and  amplitude,  intersecting  at  a  given 
ground  point.  Each  amplitude-phase  pair  is  regarded  as  a  complex 
number.  The  resultant  geometrical-optics  field  amplitude  at  such  a 
ground  point  is  given  by  the  vector  sum  of  these  quantities.  Once 
this  resultant  is  calculated,  for  both  perturbed  and  unperturbed 
profiles,  it  is  a  simple  matter  to  interpolate  the  difference  between 
both.  The  heater- induced  intensity  change,  in  dB  (relative  to 
1  V^/m^),  is  plotted  in  Figures  7a  through  7c. 

Notice  that  the  predicted  effects  occur  both  in  a  region  just 
past  the  skip  distance,  as  well  as  several  hundred  kilometers  beyond. 
Except  at  the  skip  distance,  the  predicted  field-strength  changes  are 
small,  ranging  from  about  a  dB  change  for  the  nighttime  c.ase  iFig.  7r) 
to  about  2  dB  for  the  daytime  models  (Figs.  7a  and  7c).  (Tliere  rav 
tube  cross  sections  shrink  to  zero  (that  is,  at  the  skip  distance), 
this  ray  amplitude  calculation  is  not  reliable.  Nickisch  asserts  that 
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Figure  7c.  Change  in  ground  field  intensity,  Bochkarev  model 
15  MHz,  day. 


the  power  densities  calculated  by  his  method  may  be  low  by  as  much  as 
5  dB  in  the  skip  region.  In  any  case,  the  calculated  phases  of  inter¬ 
ference  phenomena  at  the  skip  will  not  be  accurate  enough  to  support 
detailed  comparison  of  the  modified  and  unmodified.  (This  is  because 
the  path-integrated  refractive  index  is  never  a  numerically  accurate 
measure  of  plane-wa'e  phase  in  the  vicinity  of  a  caustic  [9]).  The 
second  region  of  noticeable  intensity  change  is,  however,  far  from 
both  the  interior  caustic  and  from  the  horizon  focus.  Figure  7a  shows 
a  2  dB  fade  centered  150  km  past  the  skip;  Fig.  7b  shows  about  2  dB 
750  km  past  the  skip;  and  Fig.  7c  shows  about  1  dB  500  km  past  the 
skip  distance. 

Figure  8  illustrates  tlie  steps  comprising  the  intensity-change 
calculation  for  the  daytime  curtain  antenna  model.  Notice  the  two 
amplitude  components,  corresponding  to  a  high  and  a  low  ray.  Far  from 
the  skip  distance,  there  is  only  one  such  component. 

The  remaining  figures  (9  through  12)  pertain  to  the  15  MHz 
daytime  (Bochkarev)  model  and  show  various  dependences  of  the  ground- 
level  signal. 

Figure  9a  indicates  the  trend  in  field-strength  change  as  trans¬ 
mitter  effective  power  is  increased.  As  expected,  the  greater  the 
effective  power,  the  greater  the  field  change.  Although  increasing 
power  (a  fortiori,  increasing  E^/Ep^)  much  past  the  85  dBW  level 
results  in  an  abuse  of  method  (the  differential  equations  we  solve  for 
n/NQ  are  only  valid  when  E^/Ep^  is  much  less  than  1).  one  can  infer 
that  the  90  dBW  prediction  (where  E^/Ep^  -  1)  gives  an  upper  bound  on 
effects  which  can  occur  as  a  result  of  linear  kinetic  processes 
alone.  Figure  9b  indicates  the  trend  in  this  calculation  once  we  are 
well  beyond  the  region  of  applicability  of  the  linear  theory;  here 
E2/Ep2  is  on  the  order  of  10,  and  the  predicted  ground  effect  has 
increased  by  an  order  of  magnitude. 

Figure  10  shows  the  dependence  of  the  ground-level  field-strength 
change  upon  the  local  magnetic  field  dip  angle.  The  direction  of 
propagation  has,  herein,  been  taken  as  east-west.  The  diffusion 
decreases  as  the  dip  angle  tends  toward  the  horizontal,  because  the 
gradients  in  the  field  quantities  in  Eqs.  (U)  and  (18)  become  essen- 
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Figure  9a.  Ground  effect  versus  power-gain,  Bochkarev  irodel ,  15  MHz, 
day,  60®  dio . 
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Figure  9b.  Ground  effect  versus  power-gain,  Bochkarev  model,  15  MHz, 
day,  60°  dip.  (Trend  illustrating  prediction  of  linear 
theory  in  nonlinear  regirre,  95  dbW.j 
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Figure  10.  Ground  effect  versus  dip  angle.  Bochkarev  model,  15  MHz 
day,  85  dBW. 


tialLy  perpendicular  to  the  magnetic  field  (and  are  determined  only  by 
the  azimuthal  spread  of  the  transmitter  pattern) .  The  maximum  mag¬ 
nitude  of  u/Nq  will  therefore  tend  to  increase.  However,  since  the 
vertical  dimension  of  the  greatest  perturbations  will  become  small 
with  decreasing  dip  angle,  one  must  not  assume  that  the  ground  effect 
«  will  be  monotonic  in  response.  In  fact,  some  complex  relationship 

holds  between  the  maximum  magnitude  of  n/NQ,  the  size  of  the  region 
over  which  the  density  change  diffuses,  and  the  signal  change  on  the 
ground.  This  behavior  is  evident  in  the  parametric  study  (Fig.  10). 
The  ground  effect  near  the  interior  skip  increases  with  decreasing  dip 
angle,  but  near  the  horizon  it  begins  to  decrease.  There  is  also  a 
competing  trend  in  east-west  propagation.  As  the  field  lines,  and 
hence  the  diffusion  direction,  incline  more  toward  the  horizontal, 
rays  will  tend  to  refract  at  right  angles  to  the  east-west  merid’^.i. 

In  our  calculations,  however,  we  have  not  allowed  for  ray  refraction 
out  of  the  plane  of  propagation:  we  have  taken  the  perturbed  region  to 
be  cylindrically  symmetric. 

Lastly,  Figs.  11  and  12  show  how  the  effect  of  the  perturbation 
on  ray  arrival-angle  varies  with  the  power-gain  product  and  the  local 
magnetic  dip  angle,  respectively.  In  the  first  case,  we  have  fixed 
the  dip  angle  at  60°  and  varied  the  power-gain  product;  in  the  second, 
the  power-gain  is  fixed  at  85  dBW  and  the  dip  angle  is  varied.  The 
examples  indicate  only  a  few  minutes  of  a  degree  change  in  ray  angle 
of  arrival. 
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Figure  12.  Change  in  arrival  angle  versus  dip  angle,  Bochkarev  model, 
15  MHz,  day,  85  dBW. 


IV.  CONCLUSIONS 


Our  calculations  indicate  ground-level  field-strength  changes  of 
several  dB  might  be  produced  by  joule-heating  of  the  ionosphere  bv 
intense  oblique  HF  waves.  Although  small,  these  calculated  field- 
strength  changes  are  similar  to  the  3  dB  change  measured  by  Bochkarev 
and  associates  [1,2, 3, 4),  who  gave  only  sketchy  details  of  their 
experimental  parameters.  Our  calculated  arrival-angle  change  of 
several  arc  minutes  are  much  smaller  than  the  2  to  3  deg  change 
reported  by  Bochkarev.  Our  results  are  extremely  sensitive  to  the 
model  ionosphere  used  and  therefore  indicate  that  an  experiment  should 
employ  the  widest  possible  range  of  frequencies  and  propagation 
conditions.  An  effective  power  of  90  dBW  is  far  more  likely  to 
produce  a  detectable  signal  change  than  a  power  of  85  dBW.  Our  con¬ 
clusions  are  based  solely  on  joule-heating  from  linear  waves,  and  omit 
plasma  instabilities. 

Our  results  pertain  to  steady-state  effects.  The  time  scale  for 
these  phenomena  is  roughly  the  characteristic  diffusion  time  of  the 
ionospheric  layer  in  which  strongest  electric  field  focusing  has 
occurred.  This  quantity,  taken  from  Table  2,  Vol.  1,  is  a  func¬ 

tion  of  the  ambient  recombination  coefficient  and  plasma  density 
(Eq.  [19]).  An  experiment  which  took  into  account  this  waiting  time 
might  observe  dipping  and  peaking  of  signal  strength  just  past  the 
skip  distance  (Figs.  7a  through  7c).  The  difficulty  of  such  measure¬ 
ments  is  that  the  diffusion  "waiting  time”  is,  in  some  cases,  con¬ 
siderably  longer  (10^  -  10^  s)  than  the  average  period  of  ambient 
ionospheric  fluctuations. 
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